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General stability criterion of two-dimensional inviscid parallel flow 
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General stability criterions of two-dimensional inviscid parallel flow are obtained analytically for 
the first time. First, a criterion for stability is found as u _ u > — fii everywhere in the flow, where 
U a is the velocity at inflection point, fii is eigenvalue of Poincare's problem. Second, we also prove a 
principle that the flow is stable, if and only if all the disturbances with c r = Us are neutrally stable. 
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Finally, following this principle, a criterion for instability is found as tjzjq < "Mi everywhere in 
the flow. These results extend the former theorems obtained by Rayleigh, Tollmien and Fj0rtoft 
and will lead future works to investigate the mechanism of hydrodynamic instability. 

PACS numbers: 47.15.Fe, 47.15.Ki, 47.20.Cq, 47.20.-k 



The stability of flow is one of fundamental and the 
most attracting problems in many fields, such as fluid 
dynamics, astrodynamics, oceanography, meteorology et 
al. The general stability criterion of inviscid parallel flow 
is very important for both theoretic research and applica- 
tion. A general way to study the problem is investigating 
the growth of linear perturbations by means of normal 
mode expansion, which leads to famous Rayleigh's equa- 
tion 0. Using this equation, Rayleigh first proved the 
Inflection Point Theorem, which is a necessary criterion 
for instability. Later, Fj0rtoft.3| found another neces- 
sary criterion for instability, which is a stronger version 
of Rayleigh's Theorem. These theorems are well known 
and are applied to understand the mechanism of hydro- 
dynamics instability 0, 0, Q ■ Unfortunately, both the- 
orems are necessary criterions for instability, neither of 
them can give a sufficient condition for the instability. 
And a sufficient condition for the stability is also needed 
for applications. The aim of this letter is to find a crite- 
rion for these flows. 

For this purpose, Rayleigh's equation for an inviscid 
parallel flow is employed 00013- For a parallel flow 
with mean velocity U(y), the streamfunction of the dis- 
turbance expands as series of waves (normal modes) with 
real wavenumber k and complex frequency lu = u r + iuji, 
where loi relates to the grow rate of the waves. The flow 
is unstable if and only if u>i > 0. We study the stability of 
the disturbances by investigating the growth rate of the 
waves, this method is known as normal mode method. 
The amplitude of waves, namely <f>, holds 



{U-c)(<p" -k 2 (f))-U"(j) = 0, 



(1) 



where c = u)/k = cv + ici is the complex phase speed. 
The real part of complex phase speed c r = u) r /k is the 
wave phase speed. This equation is to be solved subject 
to homogeneous boundary conditions 



= at y = a,b. 



(2) 
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It is obvious that the criterion for stability is loi = 
(cj = 0), for the complex conjugate quantities 4>* and c* 
are also the physical solution of Eq. Q . 

From Rayleigh's equation, we get the following equa- 
tion 



l^'l| 2 + fc 2 ||0l| 2 ) + 



U" 



U"(U-c r ) 



cf>\\ 2 } dy = 



(3) 



U-c 



:Hfdy = 0. 



Rayleigh and Fj0rtoft proved their theorems by Eq.©, 
which still plays an important role in the following dis- 
cussion. To find a stronger criterion, we need estimate 
the rate of J \\4>'\\ 2 dy to J ||^>|| 2 ciy. This is known as 
Poincare's problem: 



W\\ 2 dy = a 1 



2 dy, 



(4) 



where the eigenvalue /x is positive definition for 0^0. 
The smallest eigenvalue value, namely hi, can be esti- 
mated as fii > (f^) 2 - 

Then there is a criterion for stability using Poincare's 
relation (0}, a new stability criterion may be found: the 



flow is stable if 



u" 



u-u s 



> — ji\ everywhere. 



To get this criterion, we introduce an auxiliary function 
f(y) = xj-xj i where f(y) is finite at inflection point. We 
will prove the criterion by two steps. At first, we prove 
result 1: if the velocity profile is subject to f(y) > — fii, 
then cy can not be U s . 
Proof: Otherwise, 



Mi < 



U" 



U"{U-U S ) ^ U"{U -U s ) 



< 



u-u s (U- u s ) 2 - (U - Us) 2 + cj 

and if c r — U s , this yields to 

U"{U-U S ) 



(5) 



[(110' 



^ihi 2 ) 

[(Mi 



\\u-4 2 
U"{U - u 
\\U-c\\* 



dy > 
2 ] >0. 



(6) 
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This contradicts Eq.JSJ). So result 1 is proved. 

Then, we prove result 2: if — /xj < f(y) and c r ^ U s , 
there must be cf = 0. 

Proof: Otherwise if cf ^ 0, so according to Eq.J3J), for 
any arbitrary number Ut € K which does not depend on 
y, it holds 



l</>'H 2 + fc 2 |H| 2 ) + 



U"(U-Ut) , 
\\U-c\\ 2 



2 }dy = 0. (7) 



But the above Eq.JJJ) can not be hold for some special 
Ut- For example, let Ut = 2c r — U s , then there is (U — 
U s )(U-U t ) < \\U-c\\ 2 , and 

u '' (u - u it = m (n -"; )[u -, Ut) >-^ (8) 



\\U-cP 
This yields to 



{ll^'f + Ifc 2 



\\u- 



U"(U - U t ) 



\\u 



}H\\ 2 }dy>0, (9) 



which also contradicts Eq. Q . So the second result is also 
proved. 

Using 'result 1: if f(y) > — /ij then c r ^ C/ s ' and 'result 
2: if /(y) > — /Zi and c r ^ U s then c, = 0', we find a sta- 
bility criterion. Theorem 1: If the velocity profile satisfy 
> — Hi everywhere in the flow, it is stable. This 



u" 
u-u s 



criterion is more powerful than Fj0rtoft's Theorem. As 
known from Fj0rtoft's theorem, the necessary condition 
for instability is that the base vorticity £ = U' has a lo- 
cal maximal in the profile. Noting U" /{U — U s ) ~ 
near the inflection point, where £ s is the vortictiy at in- 
flection point, it means the base vorticity £ must be con- 
vex enough near the local maximum for instability. As 
shown in Fig^ there are three vorticity profiles, which 
have local maximal at y — 0. Profile 2 (U = sin(7ry/2)) is 
neutrally stable, while profile 1 (U = sin(y)) and profile 
3 (U = sin(2y)) are stable and unstable, respectively. 

Both Theorem I and FigQ]show that it is the vortic- 
ity profile rather than the velocity profile that dominates 
the stability of the flow. This means that the distribu- 
tion of vorticity dominates the shear instability in parallel 
inviscid flow, which is of essence for understanding the 
role of vorticity in fluid. In fact, we can control the hy- 
drodynamic instability just by controlling the vorticity 
distribution according these results. This is an very fas- 
cinating problem, but can not be discussed here. These 
results may shed light on the investigation of vortex dy- 
namics. 

If the flow is unstable, there are some unstable distur- 
bances with positive growth rate. For each wavenumber 
k, the fastest growing disturbance is the most important 
one. So, what the real phase velocity c r is for the fastest 
growing disturbance? Here we will show that the real 
phase velocity c r of the fastest growing disturbance is 
just the velocity at inflection point U s . 

Result 3: If the velocity profile is subject to — pp,\ < 
f(y) < — Mi with 1 < p < co, then fastest-growing distur- 
bance with highest Cj must have phase velocity c r = U s . 
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FIG. 1: vorticity profiles, profile 2 ( £ = cos(-7rj//2), dashed) 
is neutrally stable, while profile 1 (£ = cos(y), solid) and 
profile 3 (£ = cos(2y), dash doted) are stable and unstable, 
respectively. 



Proof: Suppose U s = in Eq.Q, this means Galilean 
transform for U, which make no difference for the follow- 
ing discussion. If the flow is unstable, it holds 



U(U-Ut) 1 
\\U-c\\ 2 > v 



(10) 



Otherwise the flow must be stable according to Theorem 
1. From Ea. ijTU|) . there is 



u 2 



u t )u> 



(11) 



Since Ut is arbitrary, so the left of the inequality is irre- 
spective to c r . We can rewrite the inequality as 



< P (U 2 - U t U) 



(12) 



Obviously, cf reach its largest value at c r = 0. So there 
is c r — U s , for U s is zero. 

Result 3 is very important for understanding the insta- 
bilities of inviscid flows. Since c 2 reach its largest value 
at c r = U s , the most unstable disturbances will propa- 
gate with phase speed of U s ■ If all the disturbances with 
c r = U s are neutral stable, then the other disturbances 
must be neutral stable too. This conclusion is of essence, 
and we state it as a principle. 

Principle: The flow is stable, if and only if all the dis- 
turbances with c r — U s are neutrally stable. 

By this principle, Rayleigh's criterion can be obtained 
easily. Since there is no inflection point, there is no dis- 
turbance with c r = U s according to Howard's Semicircle 
Theorem. Then the flow is stable according to the Princi- 
ple, which is what Rayleigh's Theorem states. Theorem 
1 obtained above can also derived from this Principle, 
given result 1. The criterion for instability can also de- 
rived from this Principle, we state it as a new theorem. 

Theorem 2: If the velocity profile is subject to f{y) < 
everywhere in the flow, the flow is unstable. 
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We prove this theorem by proving the following result. 
Result 4: If the velocity profile is subject to f(y) < — /i X 
everywhere in the flow, at least one of the disturbance 
with c r = U s is unstable. 

Proof: According to Eq. , for wavenumber k — and 
its first eigenfunction (f>i, it holds 

/m 2+ (TOt> » 2 ^°°- <13) 

Then cf must larger than 0. Otherwise, 

("tf^"* 1 '* (14) 

= / (A'i + Z7 ^-)ll0i|| 2 d2/<O, 

which contradicts Ea. l|13|) . So the result is proved. 

Result 4 shows that if f(y) is less than — \i\ everywhere 
in the flow, there must be unstable disturbances, so the 
mean flow is unstable for this case. This is a refinement of 
Tollmien's and Lin's who proved the similar result 
subject to the velocity profile U(y) is either symmetric 
or monotone. 

To show the power of criterions obtained above, we 
calculate the growth rate of two similar velocity profiles 
Ui = sin(1.5?/) and U 2 — sin(1.6y) with —l<y<l. 
Then there is /ii = 7r 2 /4 for estimation, and the values 
of auxiliary functions are —2.25 for U\ and —2.56 for U2, 
respectively. So U\ is stable according to Theorem 1 and 
U2 is unstable according to Theorem 2, respectively. In 
fact, there are three inflections in the velocity profile Ui- 
The growth rate of the profiles can be obtained by Cheby- 
shev spectral collocation method with 200 collocation 
points, as shown in Fig[21 It is obvious that Cj — for 
Ui and Ci > for U2, which agree well with the theorems 
obtained above. While Fj0rtoft's criterion can not point 
out the different stability between the two profiles. 

On the other hand, Arnold [lOt [Tl| discussed the hy- 
drodyanmic stability in a totally different way. He in- 
vestigated the conservation law of the inviscid flow and 
found two nonlinear stability conditions by means of vari- 
ational principle. 

Apply Arnold's First Stability Theorem to parallel 
flow, a stable criterion is < Cx < (U - U s )/U" < 
C2 < 00 everywhere in the flow, where C\ and C2 are 
constants. This corresponds to Fj0rtoft's criterion for 
linear stability, and is well knownQ- Here we find that 
Theorem 1 proved above corresponds to Arnold's Second 
Stability Theorem, i.e., a stable criterion is < C\ < 
— (U — U s )/U" < C2 < 00 everywhere in the flow. So 
It is very interesting that the linear stability criterions 
match nonlinear stability theorems very well. 

One may note that the criterions of different theo- 
rems are different from each others, U"(U — U s ) for 
Fj0rtoft's, (U-U s )/U" for Arnold's and U"/(U-U a ) for 
present works. But this make no difference for Fj0rtoft's 
and Arnold's criterions in the coarse cases, for example 
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FIG. 2: Growth rate as an function of wavenumber k, (a) 
for Ui — sin(1.5j/), (b) for U2 = sin(1.6y), both within the 
interval — 1 < y < 1. 



U"(U -U s )>0 is same as (U - U s )/U" > 0. Since the 
constants C\ and C2 in Arnold's criterions lack efficiency 
estimation for applications, they are not widely used as 
the linear criterions (eg. Rayleigh's criterion) be. 

An interesting question is which one would be the 
marginal criterion, U"(U-U S ) or U"/(U-U S )7 It can be 
seen from Eq.Q that the stability of profile U(y) is not 
only Galilean invariance of U(y), but also magnitude free 
of U(y) due to linearity. Since the value of U"(U — U s ) 
is only Galilean invariance not magnitude free, it can 
not be the the marginal criterion. While the value of 
U"/(U — U s ) satisfies both conditions, this is the reason 
why the criterions in above theorems are the functions of 
U"/(U-U S ). 

In summery, the general stability and instability crite- 
rions are obtained for inviscid parallel flow. Those results 
extend the former theorems proved by Rayleigh, Tollmien 
and Fj0rtoft. The criterions highlight the vorticity pro- 
file for understanding the instability of the parallel in- 
viscid flow. According to the criterions, the marginal of 
instability is limited in a very small zone. The criteri- 
ons for stability of parallel inviscid flow can be applied 
to barotropic geophysical flow, like Kuo did This 
extension is trivial for the cases of /-plane and /3-plane, 
and is omitted here. In general, these criterions will lead 
future works to investigate the mechanism of hydrody- 
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namic instability, and shed light on the flow control and 
investigation of the vortex dynamics. 

The work was original from author's dream of under- 
standing the mechanism of instability in the year 2000, 



when the author was a graduated student and learned 
the course of hydrodynamics stability by Prof. Yin X-Y 
at USTC. 
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